Convergence rate for numerical computation of the lattice Green's function.

Flexible boundary-condition methods couple an isolated defect to bulk through the bulk lattice Green's function. Direct computation of the lattice Green's function requires projecting out the singular subspace of uniform displacements and forces for the infinite lattice. We calculate the convergence rates for elastically isotropic and anisotropic cases for three different techniques: relative displacement, elastic Green's function correction, and discontinuity correction. The discontinuity correction has the most rapid convergence for the general case.